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ABSTRACT 


High powei density, high beam diameter, and umfoim laser beam pro- 
cessing is treated in one dimensional approximation by considering Stefan con- 
dition at the moving melt boundi\ The exact turbulent behauour of laser 
melted /one m laser cladding is not known Hence modelling is done for 
the two extreme cases complete tuibulence and no turbulence The \anation 
ol melting rate with power densitr thennal diflusmtc and clad height aie 
studied foi both models and compaied A highei powei density gnes higher 
melting rate and is rnoie economical Maximum power densitr is limited bv 
suiface melt evopoiation The turbulent melting model is more realistic than 
the stagnant melting model and hence the actual melting time is expected 
to be closer to, but somewdiat more than the melting time predicted by the 
turbulent melting model 
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Thermal diffusivity 
Specific heat capacity 
Concenterated eneigy flux 
Clad thickness 
Normalising length variable 
heat piopagation length 
Heat of fusion of clad 
Powei density 
Depth of melt 
absoibtion length 
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Nonuniformity m spatial distribution of Laser beam 

Intel action time 
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Melting temperature of clad 
Cartesian coordinate 
Density 

Normalised temperature 
Thermal conductivity 
Normalised clad height 
Normalised interaction time 
Normalised grid interval m z direction 
Normalised time interval 



Chapter 1 


INTRODUCTION 


Clad - metal composites haung a number of \aluable, parti} unique properties 
ha\o assumed ei mcioasmg importance in modem erigineeung Indmdual 
metals and aliens can not satish an\ rnoic the evei increasing stringent, com- 
plex and often c outlasting demands that materials have to meet m regaid to 
st length concision wear, and/oi heat resistance, electrical and thermal prop- 
eities etc Clad metal composites can be designed to satisfy different combi- 
nation of desired pioperties The\ can also be used as substitutes for more 
expeusn e allots The sating of scaice materials, cost reduction are among 
the other bemfits that are incurred Conventional methods of producing clad 
composites aie roll bonding, extiusion, diffusion bonding and powder metal- 
luigt routes These methods use heat and pressure to obtarn a sound rnter- 
facial diffusion bond These methods require controlled atmosphere and take 
a long tune and thus reducing the productrvity Experiments performed on 
lasei cladding by V M Weciasmghe and W M Steen [1], Y Liu et al [2, 3], 
G Coquel olio M Collm J L Fachmetti [4] W Cerri et al [5], J Com-Nougue 
and E Kerrand [6], J Singh and J Mazumder [7], RM Machmtyre [8], G 
\bbas and D R F West [9] Mmg-Chang Jeng et al [10], V M Weerasmghe, 
\\ M Steen, and D RF West [11], S J Mathews [12], E Lugscheider, et al 
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[13], Ya D Kogan [14], show that laser is a powerful tool for cladding Laser 
cladding is a process wheie mtei facial fusion bond is acheived This fusion 
bond is stiongei than the diffusion bond The time for cladding with lasei is 
much lowei than that of otliei methods and the atmosphere contiol is geneially 
not lequued 

The mam featuie of a concentrated energy flux such as laser is it’s ability 
to delivei high powei pei unit aiea (10 2 -10 4 W/mrii 2 ) to local regions on a 
woikpiece and thus lapidlv change the taiget tempeiatuie Recent advances 
in lasei technology have diawn consideiable mteiest m using lasers foi man- 
ufacturing Lasei technology not only makes manufactunng piocess snnplei 
and moie econonncal but also pioudes a unique yvay of modifying surface 
chemist ly and stiuctuie of mateiials Expenmental studies confnm that the 
CEF can In a yen useful tool m many potential applications lmohing sui- 
faco hai cloning, suiface cladding, welding cutting, suiface alloying and shock 
haidenmg Hoiuyei foi it s use to be consolidated it is necessary to pi edict 
how it yyould yyoik m noycl situation^ and noyel materials In unclei standing 
hoy it yyoiks the expenmentei is faced yyith multipaiameter problem which 
is difficult to sohe yyithout extensive factorial experimentation The principal 
vanables aie the substiate thermal and optical pioperties, the laser power, 
powei distribution and mtei action time Alternatively, an assumed physical 
pictuie of the piocess can be mathematically modeled and the model’s result 
compaied to expenmental result to prove the model’s validity and thus by 
mfeionce the physical model 

A model capable of predicting expenmental result means that previously 
unmeasui able paiameteis can be estimated In metal cladding the most im- 
portant data aie (1) the thermal cycle in each location in the fusion and heat 
affected zone which defines the extent of any phase change, (2) the peak temper 
atuie distribution and (3) the cooling lates which would affect the formation 
of metastable stiuctuies such as martensite Without a mathematical model 
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these values are very difficult to obtain 

The objective of the present woik is to model the surface cladding pro- 
cess, wlioie the clad metal powder of uniform thickness is placed over the 
substiate and is irradiated with a stationary laser source so that the clad 
melts completelv and foims a strong metallurgical bond with the substrate 
The model should predict the melt boundrj velocity, the peak temperature 
at all points and the temperature distubution Uniform, high beam diameter, 
high beam mtonsits laser cladding piocess is modeled in one dimensional ap- 
pi emulation In considering the nonlinear Stefan condition at the moung melt 
boundaiv Since the exact turbulent behauour of the melt is not known, the 
two extieme cases 1 e complete turbulence and no tuibulcnce aie considered 
The tempeiatuie chstiibution dunng pieheating and the preheating time are 
obtained liom the' Laplace tiansloim solution ol the governing one dimensional 
heat equations and boundary conditions The temperature distribution during 
melting and the melting time are obtained bj solving the one dimensional heat 
equations along with the boundan conditions and the energv balance equation 
b\ Ciank - Nicholson finite diffeience method Cladding of nickel and stainless 
steel mei C - 45 steel at clad heights of 0 2 mm 0 5 mm and 1 mm and at 
various powei densities are considered for parametric study 


Chapter 2 


LITERATURE REVIEW 


2 1 Laser Material Interaction 


\\ lieu the lasoi falls on a taiget suiface, part of the energ'v flux is reflected and 
a pait is absoibed m a thin finer The eneigj of absoibcd light is gi\en upto an 
election gas via the tuple collision of an electron, a photon and for example a 
lattice defect, which as letealed b\ calculation [15], is a necessary condition for 
fulfilling the law of conservation of momentum The collision cause electron gas 
tempeiatuie to glow, but the lattice temperature generally remains invariable 
for a definite time, because of laige difference m mass between electron and ion 
Estimates disclose that the equalisation of electron gas temperature occurs for 
about 10- 13 -10- 17 sec [15], which is much shorter than the time of energy 
tiansfei from electrons to ions 

The electron to ion energy transfer corresponds to electron photon in- 
teraction, which tend to equalize the temperature of the electron gas and the 
lattice This piocess lasts around 10 -n sec , so that the notion of the heat 
source is valid after a dwell time of 10~ 9 to 10 -1 ° sec 


When a triangular pulse interacts with a material, the maximum tem- 
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peiatuic diffeience between the electrons and the lattice is given by [15] 

(T„ - T t ) max = 

as/irar 

whe re 


T c - electron gas temperature 
T t - ion temperature 
q - flux density 

a - electron lattice heat transfer coefficient 


r - light pulse length 
a - thermal diffusivity 

At q = 10 1 H / mm 2 and r = 10 -i s, (T t — Ti) max < 1 A at q — 10 7 W/mm 2 
and r = 10" u s, (T t -T,) maa = 2000 A In both cases, a = 6 x 10 6 W / mm 3 6 
and « = 100 s In lasei beam material processing the pulse duration is 

nuuli longei than 1()~ CJ s and the powei density is lowei than 10' 11 /mm 
so tin notion of a heat souice is justifiable So the practical tasks concerned 
with CEF matenal heating at pulse duiation m excess of 1CT 9 s and at power 
densiU lowei than 10 7 W/mm 2 can m the mam be considered on the basis of 
the mechanisms of common heat conduction, be it linear or nonlinear 


The heat produced by lasei interaction is sufficient to cause change of 
phase and surface properties With very high flux or longer heating times, 
evaporation will occur at the surface, often accompanied by the ejection of 
molten matenal Still higher heating rates will produce gross vapourisation, 
the high velocity luminous plumb of vapour and particles and sometimes the 
propagation of potentially damaging stress waves These can result m the emis- 
sion from the surface of highlv energetic ions and electrons and the formation 
of plasma plume that retreats from the surface at rates as high as 10 8 mm/s 
The density of the plume ma\ be so much that subsequent laser is absorbed 
bv plasma plume rather than b\ surface Applications using continuous lasers 
or involving welding or drilling of bulk materials generally require heat fluxes 
that are well below those that cause the formation of plasma 



2 2 Uniqueness of laser source 


6 


2.2 Uniqueness of laser source 


Intel ost m lasei as a processing tool is as a xesult of it’s uniqueness as a source 
ol heat 

1 The optical eneigv is conveited into heat in a shallow legion The ab- 
soiption depth may be the mder of 20 nm 

2 The high spatial coheience of lasei beams allow focussing the beam to 
small spot of high optical intensity (10 4 - 10 8 ) W/mrn 2 This allows heat- 
ing the maternal on a \ci\ fine spatial scale bringing the temperature of 
the heated maternal m a very short time The thermal penetration depth 
is small and this piopertv combined to small optical penetration makes 
lasei heating piimanh a sm face phenomenon 

3 The lasei can be used without anv physical contact with the work piece, 
theieby minimising possible contamination 

4 Lasei beams are easih modulated electronically This allows interfacing 
the lasei into complex computer controlled processes 

The ability to effectively exploit these characteristics is enhanced by the 
availability of a number of different types of laser that are suitable for material 
processing The choice of a laser depends on the particular application and the 
material to be pi ocessed More than one type of laser will offer the combination 
of parameters necessary to do a job, and a decision to use a particular laser will 
be made on the basis of cost, a\ailabiht\ of equipment and compatibility with 
other aspects of the fabrication sequence The most commonly used lasers are 
pulsed ruby and Nd-glass, and pulsed and continuous C0 2 lasers 
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2.3 Laser Cladding 


Acheivement of a sound interfacial fusion bond is the main aim of laser cladding 
process Laser cladding offers some unique advantages 

• Conti oiled dilution le\el 

• Localised heating v Inch reduces thermal distortion 

• Controlled shape 

• Good thermal bonding 

• Fine quench microstiuctuies 

These advantages usualh lead to fundamental economic advantage ovei com- 
peting piocesses due to reduced after machining costs, and reduced product 
l ejection due to unacceptable distortion Added to these advantages are the 
improved metallurgy of the clad layer due to it’s moderately rapid quench and 
it s ease of application since the process can be a noncontact one 

2 3 1 Methods of Laser Cladding 

Laser cladding can be acheived m the following ways The simplest way is 
to preplace the metal powder over the substrate with or without the help of 
a binder and expose to a lasei source so that the powder melts completely 
and forms a strong fusion bond with the substrate (fig 2 1) The laser source 
c an be moving oi stationary In order to cover areas considerably larger than 
the diametei of the lasei beam, successive partially overlapping passes are 
deposited with a moving lasei source The other way is to inject powder 
particles by pnuematic povv der delivery system into the molten pool produced 
bv the laser beam (fig 2 2) 



2 3 Laser Cladding 


8 




Figure 2 1 Lasei cladding with preplaced powder 




Figure 2 2 Lasei cladding with pnuematic powder delivery 


2 3 2 Clad substrate combination 


Selection of a clad-substrate combination mainly depends on it’s end use Most 
metal combinations can be clad provided the clad material deoesn t have a sig- 
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Laser beam Center 



mhcanth liiglun melting point than the substiate [1] The mixed solidification 
cliiection m these cases lead to interface poiosity The buttle compound for- 
mation at the mteiface must be avoided The clad and substiate material 
should ha\e sufficient ductiliti to compensate for thermal stresses [1, 2] 

2 3 3 Experimental Results 

V M Weei asmghe and W M Steen [1] obtained very good results for cladding 
316 stainless steel on En 3 mild steel, stellite on mild steel, stellite on tool 
steel (TiC), and stainless steel on mild steel, with a moving laser beam and 
pnucmatic powder delivery svstem 'll Liu, J Mazumder and K Shibata [2] 
ueie succesful m cladding bronze over A. 4.333 A1 alloy mspite of the signifi- 
cantly highei melting temperature of bronze (1300 K) than 4.1 alloy (800 K), 
with a moving laser beam and pnuematic poyvder delivery system 4. Liu et 
al [3] obtained good results for cladding \i over AA333 A1 alloy with Cu as 
the intermediate layer, with a moving laser beam and pnuematic powder deliv- 
er} system J Singh and J Mazumder [i] obtained excellent vear properties 


2 3 Laser Cladding 


10 


for laser clad Cr, Mn, and C over AISI 1016 steel G Coquerelle et al [4] 
obtained good results for cladding Cobalt and nickel stellite on the edge of 
turbine blades W Ceiri et al [5] obtained hard, and wear resistant Fe-Cr- 
Al-Y alloy + TiC oi WC clad over AISI 310 stainless steel using a stationary 
Gaussian somce and pnuematic powdei delivery system J Com-Nougue and 
E Ken and [6] obtained good clad characteristics for cladding Co-Cr alloy over 
12% Cr steel using a 3 KW laser, moving laser beam and pnuematic powder 
delivery system J Singh and J Mazumder [7] developed m situ Ni-Cr-Al-Hf 
alloy over AISI 1016 steel using a 10 KW lasei and pnuematic powder deliv- 
er v R M Macluntue [8] was succesful in lasei hard-surfacing of turbine blade 
sluoud inteilocks using a 2 KW C'0 2 lasei moving carnage and pnuematic 
powdei deliveiv svstem G Abbas and DRF West [9] produced hard and 
wear resistant allo\ and composite clad lasers of SiC and Stellite on En 3b 
mild steel using a 2 KW CW T C0 2 lasei, mo\ing lasei beam, and pnuematic 
powder delivciy Mmg-Chang Jeng et al [10] investigated the wear properties 
and the wear mechanisms of stellite-6 powdei cladded over AISI 1020 steel 
A 2 5 KW C0 2 lasei, pnuematic powder deliverv and moving laser beam was 
used VM Weeiasmghe et al [11] cladded 316L stainless steel over mild 
steel using a 2 KW 7 C0 2 laser, moving substrate and pnuematic powder de- 
livery S J Mathews [12] investigated laser fusing of preplaced powder paste 
as viable hardfacmg method A variety of stellite powders were cladded over 
steel with a 1 2 KW 7 C0 2 laser with moving beam Ya D Kogan [14] analysed 
lasei surface alloying of iron, structural and tool steels and Al, Cu, and Ti and 
their alloys E Lugscheider et al [13] obtained good coatings of Co-Cr- W-B- 
Si and Go-Ni-Cr-B-Si hardfaemgs using a preplaced powder paste and moving 
substrate 

The effect of the process parameters, lasei power, laser beam diameter, 
traverse speed, powder feed rate and powder size and shape on the dependent 
vanables of clad shape, dilution, porosity, cracking, surface finish, substrate 
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distortion, cladding rate from experimental results are discussed below 

(a) Clad shape 

V M Weeiasmghe and W M Steen [1] observed three basic cross secton pro- 
files of single tracks as shown in hg 2 4 Profile 2 4 c is the preferred section 
lor overlapping tracks to cover an area At higher powder feed rates or lower 
powei densities profile as m hg 2 4 a is produced while at lower powdei feed 
lates oi higher powei densities profile as in fig 2 4 b is produced For a given 
ti averse speed and powder flow, there is an optimum spot diameter for max- 
imising the clad late Defoe cused beams produces profile as shown m hg 2 5 

(b) Dilution 

For a gnen specific eneigc the powdei mass flow conti ols the lev el of dilution 
[1] Dilution is independent of cladding speed for cladding 316 stainless steel 
ouu mild steel [l] Traveise speeds lower than 8 6 mm/s tend to produce large 
dilution foi cladding bron?e o\er AA333 A1 alloy [2] A higher laser power 
increased interface region and dilution [2] Thinner clad tracks lead to large 
dilution [2] Increasing the helium flow reduced the dilution but, a large helium 
flow increased poiositj [2] Increase in traverse speed increased dilution and 
decreased deposition rate m hardfacmg with a preplaced paste of Ni-Cr-B-Si 
clad [13] The hardness of the cladding was mfuenced by the degree of dilution 
and the dilution increased with increase of traverse speed [13] 

(c) Porosity 

Porosit) in laser clad maj be caused by one or a combination of the following 
cuMties between tiacks (interrun porosity, fig 2 6), solidification cavities, gas 
evolution Interrun porosity occured when two tracks of section profile similar 
to fig 2 4 a were laid side bv side with insufficient overlap [1] It could be totally 
eliminated by operating in such a way m fig 2 4 c It can also be eliminated 
ultiosomcallv vibrating the substrate during laser cladding or b\ angling 
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the powdei feed around 10° to the line of the clad track Solidification cavities 
occui at oi neai the interface between the cladding and the substrate if the 
clad has a significantly higher melting point than the substrate [2, 1] Inter- 
11111 poiosity occurs with high aspect ratio and large overlapping percentage 
[2] foi cladding bronze over AA333 A1 alloy The best overlapping clads were 
achened tor the above case with a combination of aspect ratio 3 to 3 7 and 
o\eilap peicentage 40-55 

(d) Cracking 

Clacking occuis due to tensile solidification stresses which result from differ- 
ential theimal expansion m the steep theimal gradients associated with laser 
pi crossing (fig 2 7) These tensile stresses result m cracking of nonductile clad 
lawns oi lavois with stiess laiscrs due usualh to poiosiU Thus ciacking is 
obsoiw'd on Fe/B but not on stainless steel oi stellite [1] The theimal stresses 
developed tor a cold substrate with no ielaxation is given as 

a = -~P=r P AT (2 2) 

hjQ T &S 

where 

E, - Youngs modulus of substrate 

E t - Youngs modulus of clad 

P - thermal expansion coefficient 

AT - temperature change 

Cracking of difficult systems can be eliminated by preheatng [1, 2, 3] 
The safe power input corresponding to a certain preheating temperature for 
cladding of Ni over A1 alloys is estimated by [3] 

P~P f < (850 - T h )C v (2 3) 


wheie 

P l - total power 

Pj - pow'er needed for fusion of Ni allo\ 
T h - preheat temperature 
C p - specific heat capacitv 
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Interdendritic cracking was observed m clad layer produced at high cladding 
speeds and with high power densities [11] 180° bend tests on stainless steel 
clad layers over En3 mild steel revealed only vertical cracks through clad layers 
oi none at all Theie were no examples of delaminating or interfacial cracks [1] 
On stellite/monel, monel/steel, bronze/steel clads a side blow with a chisel, 
very rarely removed the clad at the interface 

(e) Surface finish 

The degree of overlap of adjacent tracks and the section profile of each track 
govern the suiface finish of the clad layer With sufficient overlap (50 - 60 %) 
a good surface finish can be obtained (fig 2 8) The surface undulations caused 
by overlapping were typically less than 0 04 mm in contrast to 0 15 mm typical 
of arc weld oveilays foi 31GL stainless steel over mild steel [1, 11] 

(f) Homogemty of clad layer 

Heavily diluted labels showed an almost uniform composition, suggesting that 
there is considerable stirring m the pool The solidification mechanism is 
usually dendritic along the steepest thermal gradients A drastic change of 
compositon is observed in the interface region which is less than 0 1 mm [2, 1] 

(g) Substrate distortion 

Measurements on the HAZ indicate that there is minimal thermal penetration 
of the substrate and hence minimal distortion However there is some distor- 
tion due to contractural stress set up by the solidifying clad layer 

(h) Cladding rate and clad thickness 

Figs 2 9 and 2 10 show the relationship between the cladding speed, and the 
width and the height of a single track for a fixed beam diameter, laser power, 
and powder flow conditions [1] By varying the power, or the beam diameter 
it was shovm that the edge of the cladding has a constant pow r er density 
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The effect of the overlap on the minimum clad thickness is shown in fig 
2 11 From hg 2 8, 2 10, 2 11 all at constant laser power, beam diameter, and 
powder flow conditions, the following relationships were drawn [1] 


s = 

a — bw 

fig2 9 

(2 4) 

k = 

~T 

e 1 

fig2 ii 

(2 5) 

SH - 

d 

fig2 10 

(2 6) 

C = 

kwS 


(2 7) 


v here 

S - cladding speed (mm/s) 

iv - single track width (mm) 

H - single track height (mm) 

T - multiple track minimum clad thickness (mm) 

A - oveilap factoi 

C - coverage rate (mm 2 / a) 

o,b,(,d - constants 

Fioni the above equations the cladding rate as a function of speed and 
thickness can be obtained Bv differentiating this expression the speed giving 
the maximum cladding rate can be found and hence the optimum overlap 
required 

For cladding of bronze over AA333 A1 alloy, a rough estimation of the 
clad thickness is given as [2] 

net clad height = 0 64 P/ (2 8) 


or 


0 75 w t S p S w = 0 64 Pj 


v here 
Ul 

t 

S p 

s w 

0 64 
Pj 


width of the clad (mm) 
thickness of the clad (mm) 
traverse speed (mm/s) 
specific weight ( g/mm 3 ) 
powder utilisation coefficient 
powder feed rate (g/s) 


(2 9) 
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Here 0 75 x w \ t gives a rough estimate of the cross sectional area of the 
clad tiack So a higher powder feed rate will produce a thicker clad track For 
a given substiate condition and helium flow there is a minimum clad thick- 
ness foi which good clads can be made Using substrate with polished surface 
effectively leduced the minimum clad thickness The optimised processing 
paiameters for cladding Ni alloy on Cu plate and A1 substrate with bronze 
intermediate layer are summarized below 


substate 

Mirror 

Power 

I<W 

Speed 

mm/s 

Powder feed 
rate (g/s) 

Preheat 
temp (K) 

Cu 

Oscillating 

6 09 

15 5 

0 6833 

673 

Bion/e 

Oscillating 

29 

9 3 

0 45 

673 


A. good compact cobalt based allov clad over 12% Cr steel is obtained 
at tiat oiling speeds horn 5 mm/s to 11 G7 mm/s and corresponding to three 
different, distance between succesive passes, namely, 1 25, 1 62 and 2 mm The 
thickest coating acheived m this case is 0 96 mm at a 5 mm/s speed An in- 
crease in the overlapping rate lead to an increase m the clad thickness and 
\alues higher than 1 3 mm were reported [6] 

(i) Microstructure 

Laser cladding leads to extremely fine microstructure due to high cooling rate 
(fig 2 13) [2, 7] There is a HAZ in the overlapping clads, and m the HAZ 
there is no macroscopic seggregation [2] Uniform and fine undissolved Hf in 
the Ni matrix could be acheived by cladding with pnuematic powder delivery 
[7] Laser hardfacing of tuibine blades produces a higher quality hard-faced 
deposit [8] Microstructural studies showed excellent bonding of the SiC clad 
laver to stellite substrate [9] The SiC particles tend to seggregate to the up- 
per region of the clad as a result of it’s low density compared to stellite The 
lower region that contained the partially dissolved SiC particles is referred to as 
the stellite matrix of the clad Slighth more porosity is found m the composite 
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Fig 2 14 TiC particles in a stainless 
steel 31 6 clad matrix 


I 
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than in stellite The matrix of the clad showed a fairly uniform dendritic 
structure typical of stellite Large densities were found with lower traverse 
speeds due to high percentage of stellite 

The miciostructuie of Co based alloy clad over steel showed equiaxial 
stiuctuies surrounded by dendritic structures near the surface Some porosities 
and impurities were found m the overlapping area near the interface 


2.4 Mathematical Models 

The absorption of radiation and heating of the body can be treated m the 
unidimensional approximation when [16] 

r 0 > r n » max{e , L } 

where 

ro-radius of the lasei beam 

JViionumformity of the spatial distribution of the 
laser beam 
e -absorption length 
L = \fat -heat propagation length 


This eonditon is obviously not met m all cases of interaction of the laser 
with materials, but at large beam radius and at initial stages of heating, both 
by pulsed and continuous sources can be treated m umdimensional model 

To date, few models for laser cladding have been proposed A Kar and J 
Mazumdar [17] modeled m one dimensional approximation the laser cladding 
process where the substrate material is moved at constant speed while the 
clad material is poured onto it using a pnuematic powder delivery system and 
melted simultaneoush with a uniform laser beam as shown m fig 2 16 The 
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authors were concerned with the heat and mass transfer m the cladding melt 
ABCA to model the extended solid solution during laser cladding To carry out 
one dimensional heat transfer calculations they assumed that the pool extends 
upto infinity along AB Due to this assumption, cladding melt solidifies Along 
AB and this moves upwards along the direction of BC The Stefan condition is 
considered at the moving boundry The strip of cladding is assumed to be semi- 
cylmdrical m shape (fig 2 17) The freezing front is assumed to be planar due 
to small size of dendrites Stefan condition is used at the melt boundary The 
governing heat equations and boundary conditions are solved using the integral 
transform technique The cladding is assumed to melt almost instantaneously 
as soon as it is exposed the laser beam to reach a uniform temperature T 2 
This yields the following expression for T 2 , 


T 2 = T x + 


1 

a 


2 Pf 


Hr 2 c v 


L 


(2 10 ) 


where 

Ti - ambient temperature 

C p - average specific heat of cladding material 

p - laser power 

j - fraction of laser power absorbed by the cladding material 
r c - radius of the clad strip 
v - speed of the woikpiece 
L - latent heat of fusion of the clad material 


The results of the model were compared with experimental data Laser 
cladding was performed on Ni substrate with a mixture of Ni-Hf powdei of 
nominal composition 74% Ni and 26% Hf by w r eight The laser power, laser 
beam diameter, and the interaction time used were 5 KW, 3 mm and 1 44 s 
EPMA and STEM analysis showed that the concentration of Hf obtained m the 
nickel matrix is in excess of that predicted by the equilibrium phase diagram 
The experimental results show that 3 58 - 6 5 wt% of Hf in the \i matrix 
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The theoretical result predicted by the model is 7 15 wt% Hf m the Ni-Hf 
alloy This model was used to study the effect of laser power, interaction time, 
cladding thickness, cladding powder delivery rate on the composition of Hf in 
the Ni matrix The concentration of Hf was found to increase with increase of 
cladding thickness for a given speed of workpiece and laser power However, it 
would be erroneous to conclude that one can enrich the solid phase with solute 
to any concentration by increasing the cladding thickness indefinitely for a 
given lasei power and a given speed of the workpiece There is a critical value 
of the cladding thickness of which T 2 becomes equal to the melting point of the 
cladding If the cladding thickness is gi eater than the critical value then there 
will be some unmelted powder between the substrate and the cladding melt 
and hence the substiate will riot be clad It is found that more laser energ\ is 
requited to obtain an alloj of a gnen composition at a higher powxlei deh\ei\ 
rate than at a lowei pow'der delivery rate 


P H Steen et al [18] developed a model for predicting melt depth foi a 
single la-\ei material based on a two dimensional conduction balance equation 
Uniform source mo^ mg at a constant speed is used The predicted melt depth 


is given as 


where 

d 

u 

O-t, 


VJ m 

T, 


TTl 


1 OO 


Q 

Q 

K 


qd 
A 5 J r, 


d 


* P \ 


(&m - l)Q' 


(2 11 ) 


lasei beam width m the direction of beam motion 
nondimensionahsed scanning speed 
scanning speed 

thermal diffusn it\ of substrate 
scaled melting temperature 
melting point of substrate 
farheld temperature 
nondimensionahsed power density 
absorbed power density 
thermal conducts lty of the substrate 


The melt thermal properties and the latent heat of fusion are not considered 
This is valid for melting shallow' region 
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Lasei melting of a single plate with a moving Gaussian source has been 
modeled by S Kou and Y H Wang [19] considering three - dimensional convec- 
tion The equation of motion, equation of continuity, and equation of energy 
are solved by finite difference method They obtained excellent agreement 
with the experimental results pei formed on 60C3 aluminium sheets of 229 mm 
x 152 mm x 3 2 min at 1 3 KW laser power This is melting of a single plate 
and can not be treated as a cladding case 


M Picasso, et al , [20] modeled cladding with a gaussian moving source 
and pneumatic powder deliveiy system m three dimension (fig2 2) The powder 
is assumed to be predeposited on the substrate and the power absorbed bv the 
workpiece, P w is used to iemelt the clad The attenuation of the laser beam 
In the powdei pai tides and the surface reflection are considered The heat 
of fusion is neglected Simple heat equations in three dimension are sohed 
In Laplace transform method The temperature T within the material, for 
Gaussian distribution is given as 

P, 


wheie 


Jr 

T{x, y, z ) = T amb + / f(x, y, z, s)d. 

II ikr t J o 

x ari\ 2 


(2 12 ) 


1 -I- s' 


;exp 


— 

1 + s 2 


Ti 2s 2 / 


+ 


-) 


Tc 

p, 


amb 


a 

p 

c 

Vi 

k 

u 


ambient temperatuie 
absorbed laser power 

pCVj 

2k 

density of clad material 
specific heat of clad 
laser beam velocity 
thermal conductivity of clad 
laser beam radius 


The temperatuie T nothin the material, for cylindrical distribution of 
laser beam is given as 


r2Xl 1 

r (W )=T^ + sp ^ l -e^-Vririe 


(2 13) 
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wheie 


R 2 = (x — rcosd) 2 + {y — rsin8 ) 2 + z 1 


The melt boundary corresponds to the region which is at the clad melting 
point 


Experiments were carried out with the following data 


Laser power 

Injection angle of nozzle 
Powder \elocity 
Average powder radius 
Clad material 
Substrate 
k 


pc 

Tj 


1360 W 
55° 

2500 mm/s 
0 0025 mm 
Stellite -6 

mild steel (0 1% C) 
0 01 W/mm K 
0 0035 J /mm 3 K 
1573 I< 



The numerical results (fig 2 6) and experimental results are m good agree- 
ment 


The models discussed above correspond to laser cladding with pnuematic 
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powder delivery They are not valid for cladding with preplaced powder Sim- 
ple three dimensional models have the limitation of it’s inability to accomodate 
heat of fusion Finite difference methods involve much computation and are 
complex The present model is for preplaced powder technique and this ac- 
comodates heat of fusion This model uses less computational intensive finite 
difference method But, this model assumes one dimensional approximation, 
which is valid foi large beam diameter and small interaction time 



Chapter 3 


ONE DIMENSIONAL MODEL 
FOR LASER CLADDING 


3 1 Introduction 

The simplest wa^ of surface cladding with laser le , to preplace the clad pow- 
dei over the substrate and then expose to a uniform, stationary laser beam is 
modeled One dimensional heat transfer, which is \ahd for large beam diame- 
ter and high power density is considered The turbulence occunng m the melt 
during melting and hence the melt thermal conductivity is not exactly known 
Hence, modeling is done for the two extreme cases, complete turbulence and 
no turbulence The actual heat transfer process occunng during laser heat- 
ing are veiy complex with nonhnearities of thermal constants and anisotropic 
material propeities The laser source has changing spatial temporal charac- 
teristics As done by many others [17, 20, 22] the following assumptions are 
made throughout this work 

• Isotropic material 

• Stable thermal conductivity, density and specific heat 
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• No heat loss to surrounding by radiation 

• Ideal thermal contact at two phase interfaces 

• Complete absorption of the incident beam 

• Planar freezing front 

• All heat is absorbed at the surface only 

• The clad material has a unique melting point 

• Substrate is semi-infmite 

The anisotiopic nature of the materials depend on a number of factors 
like, pnoi treatment, and it’s spatial distribution can not be predicted The 
thcimal constants may vary to a considerable extent, but difficult to model 
Owing to a small interaction time and a large heat flux coming m the heat 
loss to the surrounding can be neglected Thermal resistance between powder 
clad and substiate is not exactly known The power density considered for 
calculation is the actual power absorbed by the clad 

3.2 Turbulent melting 

In this model the melt pool is assumed to be so turbulent that it has infinite 
thermal conductivity Consider a uniform laser source heating the surface of 
a semi-mfinite body over which the clad powder is preplaced to a uniform 
thickness The use m surface temperature is limited b\ the heat conduction 
into the material from hotter to colder metal and heat of fusion When the 
surface is below the melting point the surface condition becomes 
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When the surface reach the melting point, liquid metal forms and the melt 
boundary moves with certain velocity f into the clad Now, according to this 
model, the lasei effectively falls on the moving melt boundry One condition 


q 



Figure 3 1 Preheating with uniform source 


at the moving boundry s(t) is that 

T\ = T m (3 2) 

The othei condition at the moung boundrv is the stefau condition, i e , 

, dTi r ds 

= —p\L m — + qo (3 3) 

The one dimensional unsteady state heat conduction equation for the clad 
la\er and the substrate are gnen below 
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Figuie 3 2 Tuibulent melting with uniform source 
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Initial temperature of the clad and substrate is same througout, (T 0 ) 
At t = 0 


Ti = T 2 = To 


(3 9) 


Foi numerical calculations instead of equation(3 3) a more convenient 

foim can be obtained b\ considering the rate of change of heat m the material 

d \ r h /-co r h ft tv 

/ Pi c iTi(z, t)dz + / p 2 C 2 T 2 (z,t)dz = — P\C\T m ~ + pitj / — — dz 

(it Jo hi dt Jsu) at 


+P2C2 


L 


00 <9T 2 


dt 


dz 


d& 


— — Pi(-\T m — + X] 
dt 


dTj 

dz 


+A 2 


8T 2 


dz 


ds dds 

- PtCiTm dt p\ L m gdt 

+q{t) 


(3 10) 


Integrating the above equation over any time interval {ti,t 2 ) gives 

ft 2 rh rh 

rh . r h _ . 

+ 


/ cy(f)df = pi Ci l Ti(zM)dz - picx f Ti{z,t x )dz + 
Jt\ Js(t2) Js(Lli) 

rh rh 

p 2 c 2 / T 2 {z,h)dz - p 2 c 2 / T 2 {z,t x )dz 
JsUi) Js(ti) 


S{t2) ' Js(ti) 

Pi (ciT m + L m )(s(t 2 ) — s(ti)) (3 11) 


When the surface is below the melting point, the equations and linear 
boundary conditions can be solved bv Laplace Transform method to get the 
temperature distribution The clad laver temperature distribution is given as, 


Xl 


g n+x ierfc 

Ln =0 


2 h(n + 1) — x 
2 -J o. x t 


+ g n ierfc 

n=0 


2 nh + x 


2 y/aXpt 
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The subst rate temperature distribution is given by, 


2 a CO 

T 2 = ~ (1 + s) Yj 9 n ierfc 

A 1 n=0 


(2n + 1 )h + (x - h 
2y/aYt 


(3 13) 


where it? fc is the integral of the complementary function 


erf(x) = 
e? fc(r) — 
icrfc(r) = 


TJ^~ v dy 

1 - erf{x) 

IT erfc(y)dy 
i 

7^ 


1 e x ~ xerfc{x) 


When melting occurs, the formulated equations with nonlinear boundry 
condition can bo solved by Crank-Nicholson implicit finite difference method 
It is convenient to intioduce the following normalised \anables 


0 1 - 

T\ - To 
(T m ~ To) 

(3 14) 

9 2 = 

T 2 -Tq 
(Tm - To) 

(3 15) 

l = 

— (c,(T ra -r„) + i m ) 

Qo 

(3 16) 

c - 

2 

7 

(3 17) 

4’ = 

h 

l 

s 

1 

(318) 

£ = 

(3 19) 

T = 

a\t 

T 

(3 20) 


Now, the boundiy conditions and equations [3 2, 3 4 -3 9, 3 11 ] after normal- 
ising using equations 3 14 —• 3 20 are as below 
At ( = e(t) 


At e(t) < C < V ; 


6i = l 

86 ! 8% 


(3 21) 


(3 22) 


dr d ( 2 
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At C = 


Q\ — 9 2 
dO\ d6 2 

ac = sc 

A.t ij> < C < oo 

aj d9[ _ S 2 #i 
a 2 dr d( 2 

\t ( = oo 

0 2 = 0 


(3 23) 
(3 24) 


(3 25) 


(3 26) 


Hie eneigt balance equation becomes 

f — dr = PKi(T w -T 0 )/ T 0,(C,r 2 )dC- 

! r, a I Je{Ti) 

PiC\(T m - Tq)1 P 0i(C,r 1 )dC + 

Je(n) 
r oo 

p 2 c 2 {T m -T 0 )l J 9 2 (C,T 2 )d(- 

/•OO 

PiC 2 {T m -T 0 )l ^ ^(C.nJdC- 

p,i(ci(T m - To) + L m ){e{r 2 ) - e(n)) (3 27) 


The noimahscd equations written in the Crank_Nicholson finite difference 
foim are given below 
At c = ^(0 

fl, = 1 (3 28 ) 


At £•(/) < ( < V’ 

i9j i+i ~ - 20, n+i + ^"-i 1 + ^ n 4-i - 2gr + flr~i 


Ar 


^t C = 0 


Ai[ 


2AC 2 


0 „ + i_ c+ i _ % 


AC 


-] = A 2 [- 


AC 


(3 29) 


(3 30) 
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Figure 3 3 Energy balance during laser heating 


At 0 < C, < co 

Ai oy 1 - gr ^i 1 - w 1 + Ci 1 + ffVr - w + flu 

A r 2AC 2 


At £ — oo 


A 

imax 

( £ = o°) 

(3 31) 

(3 32) 


0 2 = O 
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The energy balance equation becomes 


oi ' 



ih — 2 

PiCi(T m - T 0 ) AC 

E «<T 

_t=istari+2 


+ 


an+l , an+l 

^(tatart+1) 


ih—2 

, E «& 

[i=ts tart-hl 


+ 


Pi c i(T m — To) AC 

ft c ;. + ^ c y Tni _ To )Ac«iU, + 


+ 


Qfistart) 0(th- 1) 


P2C2{T m — To)A( 
Pit 2 (Tj/i - To) AC 


imax — 1 

e «?,r 

i=i/i-f2 

imai-1 

E «?.) 

i=i/i+2 


+ 


fl r,An + ggu 


+ 


^ft/H-l) + ^mns) 


(3 33) 


The above eneigy balance equation is used as the basis of an iteiative 
loop Throughout the finite difference scheme a fixed value of AC is used 
While the clad is melting, the time step is taken as the time taken for the 
phase boundry to move through a distance AC Hence in one time step the 
phase boundry moves from one grid point to next 

When the surface is moving, an initial value of At is taken and the 
equations 3 28 to 3 32 are solved A finite difference approximation using 
equation 3 33 then yields a revised estimate for the time step At, this iteration 
is repeated until the desired accuracy is achieved before continuing to the next 
time step 


3.3 Stagnant melting 


This physical model assumes a stagnant melt pool of specific thermal con- 
ductivity Effectively, it will become a three layer problem Here also the 
temperature distribution during preheating is calculated using equations 3 12 
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md 3 13 as discussed in the previous section 


q 


0 



Figme 3 4 Stagnant melting with uniform source 


Heie the boundary conditions during melting become, 
At /=(), 


Ai 


&Ti 

dz 


Qo 


(3 34) 


At 0 < z < a(0» 


At 2 — s(f), 


l_dT\ = &T\ 
a\ dt dz 2 


T\ 

dTx 
1 dz 




m 

dz 


(3 35) 


(3 36) 
(3 37) 


1 dT-j _ d' 2 T 2 (3 38) 

a 2 dt dz 1 


At i>(t) < z < h, 
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At h < z < oo, 


At i = oo, 
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t 2 = r 3 
x or 2 5 t 3 

Ai 57 = A3 aT 

(3 39) 

(3 40) 

1 - d2T i 

a 3 5/ 5z 2 

(3 41) 

li 

E^ 

(3 42) 


Now, the ciK'igv balance equation becomes, 



/> 


P\L m db(i) + 


4 i 


-s(0) 


p l CiTi{z > t 2 )dz ~ 
p\C\T\(z, t\)dz + 


f ■ 

rh 


li(tl) 


p 2 c 2 T 2 {z,t 2 )dz - 
P2C 2 T 2 (^,ii)dz + 


/ Pi^Tz{z,t 2 )dz 

Jh 

r 00 

/ p^3,Ti{z,ti)dz 

Jh 


(3 43) 


The above equations [3 34 3 43] m the finite difference form, after normalising 
using equations 3 14 3 20] aie given below 

At ( = 0, 


\ /nr nr 1 aw+1 fl n +l 
Ai(i m ~ JoifVi ~ v i + 1 i 

I [ AC ' 

At 0 < C < ^(f), 

0J‘ H - 0» _ - 2ff 1+1 + + ff+l ~ ± ^-l 

„ 2AC*~~ 


(3 44) 


(3 45) 
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Ate = s(t), 


0i=l 


(3 46) 


At b(t) < C < h, 

a x 6>; H 1 - 9\ l _ 6[ 1 ^ - 2 9\ l+1 + e/ + 0,+i - 

n, At “ 2AC 2 


(3 47) 


At C — h, 


^2 ( 



) = A 3 ( 


Cl 1 - ff + * v 
At 1 


(3 48) 


At h < C < oo, 

ai fl" +I - B\ l _ 0\'lx - 2 0l l+l + 61 jj 1 + fl, w + i - 20” + 9^i 
a , At 2AC 2 ( 3 49) 


At C = 00 , 


0 = 0 


(3 50) 


The above equations [3 44 - 3 50] can be solved by the same iterative 
method used m the pievious section The finite difference method is vali- 
dated by companson with the analytical solution (equations 3 12 and 3 13) for 
pi cheat mg 


Chapter 4 


RESULTS AND DISCUSSIONS 


4.1 Introduction 


Calculations aio nude for cladding of nickel and stainless steel over C-45 steel 
The piopeities of the materials used aie listed below m table 4 1 The density 
and the specific heat of the liquid metal is not significantly different from that of 
the solid Therefore, the liquid clad layer is assumed to have the same density 
and specific heat as listed above for the corresponding solid [20], except the 
thermal conductivity The liquid thermal conductivity is not exactly known 


Table 4 1 Thermal properties of substrate and clad materials 


Properties 

0 45 %C Steel 

Nickel 

Stainless steel 

Theimal conductivity(J /s/I</mm) 
Density (g (mm 3 ) 

Specific heat capacity (J/K/g) 
Latent heat of fusion (I/mm 3 ) 
Melting pomt(K) 

Boiling pomt(K) 

Thermal diffusivity(mm 2 /s) 

0 0589 

7 849 x 10~ 3 

0 461 

2 1 
1673 
3273 

0 001628 

0 09 

8 906 x 10" 3 

0 444 

2 000 
1673 
3173 

0 002276 

0 0163 

7 317 xl0~ 3 

0 461 

2 134 
1723 
3273 

0 000483 
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The molt theimal conductivity for liquid Zn and Sn are about 50% of the 
theimal conductivity of the conespondmg solids [23] Therefore, thermal 
conductivities of liquid Ni and stainless steel are assumed to be half that of 
the conespondmg solids Initial temperature of the clad and substrate are 
considered as the loom temperature (300 K) The clad height is varied from 
0 2 mm to 1 mm 

4.2 Preheating 

The tempeiatme distribution during picheatmg is calculated using equations 
3 12 and d Id The pi cheat mg tune (till the surface reaches melting point) is 
calc ulatc'd it datively bv choosing an initial value of 0 1 s and calculating the 
clad tempeiatme at the smlace The eiror in calculated preheat time is kept 
below It) 0 s 

I lg 1 1 slums the tempeiatme distribution obtained for cladding of nickel 

ovoi steel at clad height of 1 mm and power densities of 100, 150 and 200 
M /mm 2 The tempeiatme gradient along the depth is very high The 
tempeiatme dec leases exponentially from the surface to the interior With 
me lease ol power density the temperature gradient increases and the heat 
penetration depth reduces The heat penetration is * 3 mm at a power density 
of 100 I Y/mm 2 and it drops to * 1 5 mm when the power density is doubled 
(2()0 W/mni ? ) At the clad-substrate interface a discontinuity in change of 
slope is observed because of the change in thermal diffusivity at the interface 
The tempeiatme gradient in the clad and the substrate are different Fig 4 2 
shows the temperature distribution for cladding of nickel over steel as well as 
cladding of stainless steel over steel for 1 mm clad thickness and 100 W/mm 2 
power density The temperature gradient is much higher for stainless steel 
clad due to it’s lower thermal diffusivity Hence, the heat penetration depth 
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Figure 4 1 Temperature distribution at start of melting of Ni clad over steel 
at 1 mm dad thickness 

is much smaller loi stainless steel clad (rs 0 5 mm) as compared to that for 
nickel clad ( & 3 mm) 

Fig 4 3 shows the pieheat time vs power denstitv for Ni dad at 0 2 
mm, 0 5 mm, 1 0 mm thicknesses The preheat time increases with increase 
m clad height and with decrease of power density The increase of Ni clad 
thickness me leases heat dissipation to the interior due to it’s higher thermal 
chlfusivity than steel substiate and hence the surface reaches the melting point 
at a latei time Increase of power density increases the heat input rate and 
i educes heat penetration to the interior, hence the surface reaches the melting 
pomt earlier Fig 4 4 shows the preheat time vs power density for Ni and 
stainless steel clads with 1 mm thickness and at 100 W/mrri 2 power density 
The preheating time for nickel is more than that for stainless steel clad for 
a given clad thickness and power densit\ This is due to the smaller thermal 
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Fignio 4 2 Temp (nature distribution at start of melting 1 mm thick Ni and 
stainless steel dads at 100 W/mm 2 power density 

difiusivitv of stainless steel (0 0005 mm 2 /s) than nickel (0 0023 mm 2 /s) The 
smaller thermal diffusivity leads to less heat dissipation to the interior and 
hence most of the heating occurs at the surface and thus reduces the preheat 
time The* difleience ui pieheatmg time reduces as power density increases a 
due to the sm alien theimal penetration at a higher power density 


4.3 Turbulent Melting 

Hoio, the melt pool ove. the solid clad material is assumed to be so turbulent 
that the melt pool has infinite thermal conductivity Under these conditions, 
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10 12 14 16 18 20 

POWER DENSITY ( X 10 w/mm 2 ) 

Figuic 4 3 Variation of Preheat time with power density of Ni clad at various 
clad thicknesses 

temperature of the liquid metal is maintained at the melting point For this 
case, equations 3 28 to 3 33 are solved by finite difference method as discussed 
m section 3 2 with A z as 0 004 mm The grid size independency is tested with 
A z as 0 003 mm and the difference m calculated time for the melt boundary 
to move to the next grid point is less than 10 -6 s 

Fig 4 5 shows the melt boundry migration of 1 mm thick Ni clad layer 
over steel at various power densities as a function of time It shows that the 
late of melt boundary migration increases as melting proceeds This is because 
the steel substrate has lower thermal diffusivity than the nickel clad and hence 
less heat will be dissipated to the interior substrate as melting proceeds Fig 
4 6 show's the melt boundry migration of 1 mm thick stainless steel clad layer 
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h'igmo 4 1 Variation of Pi cheat time with power density for Ni and stainless 
stool dads over stool at .1 rum clad thickness 

ovoi sit'd al vanous powei densities as a function of time Contrary to Ni 
dad, the rate of melt boundary migration for stainless steel clad decreases 
as molting proceeds This is due to higher thermal diffusivity of the steel 
substiato than that of stainless steel The higher thermal diffusivity increases 
the heat dissipation into the substrate Both figures 4 5 and 4 6 show that the 
lalo of molting ma eases with increase of power density This is due to the 
smaller heat penetration into the substiate at higher power density 

Tables 4 2 and 4 3 shows the average rate of melt boundry migration, 
the preheat time, and melting time for Ni and stainless steel clad over steel at 
various powei densities and clad thicknesses As seen from these tables, the 
melting tune is of the order of fraction of a second The tables show that the 
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lable 4 2 Results obtained for cladding of Ni over steel by turbulent model 


Power density 

'Clad 

Pi cheat 

Melting 

Average melting 

Heat input 


height 

tune 

tune 

rate 


(W/niir) 

(nun) 

00 

00 

(mm/s) 

Jjmm 3 

“Too 

""] 0 ” 

0 05(1 

0 120 

14 3 

12 0 

150 

1 0 

0 023 

0 075 

19 2 

11 3 

200 

1 0 

0 013 

0 052 

25 6 

10 4 

100 

0 5 

0 044 

0 086 

11 9 

17 2 

150 

0 5 

0 021 

0 048 

18 5 

14 4 

200 

0 5 

0 012 

0 032 

25 0 

12 8 

100 

0 2 

0 038 

0 057 

10 5 

28 5 

150 

0 2 

0 018 

0 030 

16 7 

22 5 

>00 

0 2 

0 011 

0 019 

25 0 

19 0 


Table 1 i Results obtamocl for cladding of stainless steel over steel by turbu 


a . a 

1 i ■ i \ 

i 

1 in <l 

’■ , > 

> ■ n c p • 1 . 

i i t 

(W /mm*) 

height 

(mm) 

time 

oo 

time 

(s) 

rate 

(mm/s) 

(J jmm z ) 

100 

To 

BH 

0 112 

97 

11 2 


1 0 

■ EM 

0 067 

15 9 

10 1 

fli 

1 0 

0 002 

0 050 

20 8 

10 0 

100 

0 5 

0 009 

0 075 

76 

15 0 

150 

0 5 

0 004 

0 041 

12 3 

12 3 

200 

0 5 

0 002 

0 029 

18 5 

11 6 

too 

0 2 

0 011 

0 052 

49 

26 0 

150 

0 2 


0 026 

91 

19 5 

200 

0 2 

0 002 

0 017 

13 3 

17 0 
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Pigmc* 4 5 Melt boundiv migration with time for 1 mm thick Ni clad over 
stool at various powei densities 

pi cheat time foi Ni clad is significant The preheat time for Ni clad over steel 
is 30-50 % of the total melting time whereas that for stainless steel is only 
10-15 % The total melting time for stainless steel clad over steel is always less 
than that of nickel over steel for a given clad thickness due to lower thermal 
diflusivity of stainless steel As seen from the tables 4 2 and 4 3, for a given 
powei density, stainless steel clads start melting at the same time irrespective 
of the clad thickness This is because, the heat penetration depth is less than 


the clad thicknesses 
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Piguio 1 (i Melt bouncily migration with tunc for 1 mm thick stainless steel 
dad ovoi stool at vaiiuus powei densities 

4 4 Stagnant Melting 

IIoio, the dad molt layoi is assumed to be stagnant i e no convection in the 
molt pool, and tho boat Liansfer from the melt to the clad-melt interface is due 
to only thermal conduction Foi this case, finite difference equations 3 44 to 
3 50 ale solved as discussed in section 3 2 with A z as 0 004 mm The grid size 
independency is tested with Hz as 0 003 mm and the difference m calculated 
time foi the 1 melt boundary to move to the next grid point is less than 10" 6 
s Chile illations at high powei densities predict surface temperatures much 
highoi than that the boiling point of the liquid clad layer Hence, calculations 
aie made for cladding at power densities lower than those used m turbulent 
melting model 

Pig 4 7 shows the melt boundary migration along the depth as a function 
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TIME (s) 


Figure 1 7 Melt boundry migration with time foi 1 mm thick Ni clad over 
steel at vanous power densities 

of time at powei densities of 50, 60 and 75 Wfmm 2 for Ni clad of 1 mm 
thickness ovoi steel The rate of melt boundary migration decreases as melting 
proceeds A small increase in melt boundary migration rate is obtained m the 
begming as in the turbulent melting case, but as melting proceeds the thermal 
bamei due to the liquid melt dominates and hence migration rate decreases 
Fig 1 8 shows the rate of melt boundary migration of stainless steel clad over 
steel at power densities of 10, 12, and 20 W/mm ? Melting of stainless steel 
clad staits earlier than nickel and the decrease of melting rate is more for 
stainless steel clad layer and completion of clad melting occurs later than that 
of nickel clad The time for completion of melting decreases continuously as 
power density increases The melt boundary migration rate increases with 
increase of powei density and with decrease of clad thickness Figs 4 9 and 
4 10 show the temperature distribution at the end of cladding The surface 




MELT DEPTH (mm) 



TIME (s) 

Figuio 1 8 Molt bound) v initiation with time for 1 mm thick stainless steel 
dad o\('i si ('('1 at vaiicms powoi densities 

ic'Dipo] til ui(' im leases with me nvi.se of power density 

Tables l 1 and 4 0 show the average rate of melt boundry migration, 
the p. cheat lime, moll, ms time and surface temperature for Ni and stainless 
slop] ( lulls ovol kiwi at vanous power densities The melting rate increases 
with meieasc of power density The heat input required to melt a given clad 
thickness dnei eases with meiease of power density This is because of less 
heal penel ration into the substrate at a higher power density Hence, a power 
density as high as possible should be used But, both the tables 4 4 and 
1 B show that the suifat e tcmpeiature increases with increase of temperature 
When the smfnio lemperatuio is above the boiling point of the liquid melt, 
evopoi ation o. < u.s at the surface leading to loss of metal and the cladding will 
not lie off el tive So the maximum power density that can be used is limited 



lagnant Melting 


49 


> 4 1 Results obtained for cladding of Ni over steel by stagnant melting 
4 


ioi density 

Clad 

Preheat 

Melting 

Average melting 

Heat 

Surface 


height 

time 

time 

rate 

input 

temperature 

N /mm 2 ) 

(mm) 

to 

(s) 

(mm/s) 

J/mm? 

(K) 

50 

1 

0 176 

0 769 

1 67 

38 5 

2695 

60 

1 

0 126 

0 601 

2 11 

36 1 

2883 

75 

1 

0 085 

0 449 

2 75 

33 7 

3156 

100 

1 

0 050 

0 252 

4 95 

25 2 

3561 

75 

0 5 

0 074 

0 115 

12 2 

43 1 

2694 
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:ing model 


wo i density 

Clad 

Preheat 

Melting 

Average melting 

Heat 

Suiface 


height 

time 

time 

rate 

input 

temerature 

W /mm 2 ) 

(mm) 

to 

to 

(mm/s) 

(J /mm 3 ) 

(K) 

10 

1 

1 762 1 

7 01 

0 1911 

70 1 

2872 

12 

1 

1 080 

5 71 

0 266 

68 5 

3095 

20 

1 

0 273 

3 00 

0 367 

60 0 

3979 
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Figure 4.9: Temperature distribution at the end of cladding for 1 mm thick Ni 
clad over steel at various power densities 

by surface evoporation. Hence, from the tables 4.4 and 4.5 the maximum 
power density that can be used without surface evoporation is 75 W/mm 2 
for Ni clad over steel and 12 W/mm 2 for stainless steel over steel, at 1 mm 
clad thickness. A clad melt layer with lower thermal diffusivity shows higher 
surface temperature. This model predicts a lower melting rate for a lower 
thermal diffusivity clad. 


4.5 Comparison of the Two Models 


The stagnant melting model always predict a higher melting time than the 
turbulent melting model. For a given power density and clad thickness the 
turbulent model predicts a lower melting time for a clad with lower thermal 
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Figure 4.10: Temperature distribution at the end of cladding for 1 mm thick 
stainless steel clad over steel at various power densities 

difi'usivity, For both cases, the melting rate increases with increase of power 
density. The total heat input required also reduces with increase of power 
density, since the heat penetration into the substrate is less at a higher power 
density. The turbulent melting model has no limit on the maximum power den- 
sity that can be used. But, the stagnant melting model predicts the maximum 
power density above which surface evoporation of the clad melt occurs. 

Fig 4.11 shows an example of the melt boundary depth as a function 
of time for both models for cladding of 1 mm thick Ni clad over steel at 100 
W/mrri 2 power density. It shows that melting starts at the same time. As 
melting proceeds the time gap for melting a given depth goes on increasing. 
Hence, for small melt depth the difference in calculated melt time for the 
two models is small, but for large depth the difference is significant. This is 
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Figure 4.10: Temperature distribution at the end of cladding for 1 mm thick 
stainless steel clad over steel at various power densities 

diffusivity. For both cases, the melting rate increases with increase of power 
density. The total heat input required also reduces with increase of power 
density, since the heat penetration into the substrate is less at a higher power 
density. The turbulent melting model has no limit on the maximum power den- 
sity that can be used. But, the stagnant melting model predicts the maximum 
power density above which surface evoporation of the clad melt occurs. 

Fig 4.11 shows an example of the melt boundary depth as a function 
of time for both models for cladding of 1 mm thick Ni clad over steel at 100 
W/rnm 1 power density. It shows that melting starts at the same time. As 
melting proceeds the time gap for melting a given depth goes on increasing. 
Hence, for small melt depth the difference in calculated melt time for the 
two models is small, but for large depth the difference is significant. This is 
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Figme 4 11 Comparison of melting rate of the two models of 1 mm thick Ni 
dad ovoi steel at 100 W/mm 2 power density 

because, at small melt depths the heat transfer occurmg m the clad melt is 
not smgihcaiit 

The ac tual situations would in general be somewhere in-between the ex- 
tieme conditions considered in the two models However, it is not easy to 
model the level of turbulence in the liquid and it’s effective thermal diffusivity 
The tuibulenc e occurmg m the melt would in general depend on the power 
density, the thickness of the melted layer and thermal gradients hence, it will 
also be a function of time, as melting proceeds Turbulence at the start of 
me] l mg is negligible As melting proceeds, the melt clad thickness increases 
and hence the turbulence is expected to increase As turbulence increases, ef- 
fective' thermal diffusivity of the liquid melt increases Now, more heat reaches 
the hquid/sohd interface and hence the melting occurs at a faster rate 
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The actual melting time m general would be somewhere m-between the 
melting times predicted by the two models At small clad heights and small 
power densities the melting time predicted by stagnant melting model is ex- 
pec Ic'd to lie closer to the actual melting time The difference m the melting 
time pi echcted by cither model is small at small clad thicknesses Hence, the 
melting time predicted by either model may be used as a first appioximation 
At huge dad thickness and large power densities the actual melting time is 
expected to be c lose x to, but somewhat more than the melting time predicted 
bv turbulent melting model 

Finally, at very high pciwei densities, large amount of energy is absorbed 
m a voiy thin lavei At power densities higher than a certain value, direct 
c>\ oiioi at ion may oum with little or no melting At these power densities, 
(lie dadding oi smlate moiling is not possible and even the turbulent melting 
model bleaks down 

In summary, it may be c oncluded that the turbulent melting model can be 
used to pi echc t the melting times as a first approximation, during laser cladding 
at niter mediate power densities The actual melting time is expected to be 
somewhat laigei than the melting time predicted by this model The stagnant 
melting model seems to be umeahtic, particularly at high clad thicknesses 


Chapter 5 


CONCLUSION 


Lasoi c lddchnp, plot ess with pi ('placed powdei and with uniform, stationary 
soul t o is niodelod m one dinieiisioii Since melt turbulence is not known, 
tlie two t'\t ionic’ cases, with complete turbulence and without turbulence, are 
modeled Calculations aie made for nickel and stainless steel clads over 0 45 
% t ai hem stool at vanous powoi densities and clad thicknesses for parametric 
st tide The follow nip, conc lusions are denvcd from the predicted results 

J Both models pi edict a very small melting time of the order of a fraction 
ol a second to lew seconds 

2 The pi cheat time is significant m comparison to the total melting time 
f oi clad in At dials with highei theimal diffusivity due to higher heat 
poneti at, ion into the liitenoi The preheat time for Ni clad over steel is 
'10-50 % ol the total melting tune, but the preheat time for stainless steel 
( lad is only 10-15 % of the total melting time Stainless steel has lower 

l hoi mal chllusivity 


3 Most of the heat absorbed is used for melting the clad 
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4 The heal penetration depth is m the order of a few millimeters So, the 
mtenor metal properties and microstructure will remain unaffected 

5 The total power input required to melt unit volume of the clad material 
dec teases and the melting rate increases with increase of power density 

0 The at tual melting time would in general be higher than the melting time 
pi edit ted by turbulent melting model and smaller than the melting time 
pi edicted by stagnant melting model The difference in the predicted 
melting time for both models is small at small clad thickness Hence, 
the melting time pi edicted by cither model may be used at small clad 
thickness At large clad thickness and large power densities the actual 
melting time is is expet led to be closei to, but somewhat more than the 
melting time pi edit ted bv turbulent melting model 

7 It is suggested that the tuibulent melting model may be used to esti- 
mate the melting tunes foi the clad materials during laser cladding with 
pi (‘placed powder and with uniform, stationar\ laser souice The ac- 
tual melting time is expected to be somewhat more than the predicted 
melting tune 
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